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Abstract: Based on J/i/) events collected with the BESIII detector, with corresponding Monte Carlo samples, the 
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1 Introduction 

The charged track reconstruction (tracking) efficiency 
and its systematic uncertainty play an important role 
in high energy physics experiments. A detailed study 
to develop a robust methodology for determining the 
tracking efficiency is essential for both software devel¬ 
opment and experimental data analysis. However, the 
tracking efficiency is not always given a good definition 
or straightforward determination. This paper, focusing 
on the methodology, presents a detailed study of track¬ 
ing efficiency and its systematic uncertainty, with various 
influencing factors on the efficiency, as well as a number 
of essential validations. 

The key tracking sub-detector in the BESIII detector 
[l| at the Beijing Electron-Positron Collider (BEPCII) Q 
is a small-celled, helium-based main drift chamber 
(MDC) with 43 layers of wires, which has a geometri¬ 
cal acceptance of 93% of 47r and provides momentum 
and dE/dx measurements of charged particles. The 
other components of BESIII mainly include an electro¬ 
magnetic calorimeter (EMC) made of CsI(Tl) crystals, a 
plastic scintillator tiine-of-ffight system (TOE), a super¬ 
conducting solenoid magnet, and a muon chamber sys¬ 
tem (MUC) made of Resistive Plate Chambers (RPCs). 

A clean sample of J/'i/; pp7r+7r“ decays is cho¬ 
sen for this study, due to low background level, high 


statistics, broad momentum and angular distributions, 
etc. A data sample of (255.3 ±2.^ x 10®J/-^ events col¬ 
lected with BESIII is used here. [3] The optimization of 
the event selection, the investigation of possible back¬ 
grounds, and many validations in this study are per¬ 
formed using Monte Carlo (MC) simulated data sam¬ 
ples. The GEANT4-based simulation software BOOST Q 
includes the geometric and material description of the 
BESIII detector and the detector response and digiti¬ 
zation models, as well as the tracking of the detector 
running conditions and performance. The production of 
the J/ij) resonance is simulated by the MC event gener¬ 
ator KKMC ; the known decay modes are generated by 
EVTGEN 0 with branching ratios set at Particle Data 
Group (PDG) valuesQ, while the remaining unknown 
decay modes are modeled by LUNDCHARM 

2 Definition 

The tracking efficiency {etrk) for a given particle is 
defined as 


where the denominator {N) is the number of signal 
events with all other particles in the final states required 
to be reconstructed except the one under study, no mat¬ 
ter whether the studied particle is reconstructed or not; 
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Fig. 1. Typical Pt distribution of (a) protons and (b) pions from exclusive MC sample. 


the numerator in) is the number of events with all par¬ 
ticles in the final states reconstructed, including the one 
under study. 

For example, in the calculation of the tracking effi¬ 
ciency for protons with the control sample J/4’ —t 
ppTT+TT”, at least three charged tracks should be recon¬ 
structed and identified as p, 7r+, and 7r“, respectively. 
The number of signal events {N) is obtained by fit¬ 
ting the missing mass distribution of p, 7r+, and 7r“. 

The missing mass is defined as 

where the missing four-momentum is 

determined from the difference between the net four- 
momentum of the J/ip particle and the sum of the four- 
momenta of p,7r+, and n~ in the event. The numerator 
(n) is obtained by fitting the missing mass distribution 
with the fourth charged track reconstructed. 

It is an advantage to obtain both n and N by fitting 
the missing mass distribution, since one can use the same 
function to describe the signal and part of the errors in 
the fitting can be cancelled. 

The systematic uncertainty of tracking is 

defined as the difference of tracking efficiency between 
MC {etrk{MC)) and data {etrk{data)) 

Two hundred million inclusive J/ip MC events were 
used to investigate the possible backgrounds from J/ip 
decays. Table[T]shows a topology analysis of the inclusive 
MC sample after event selection which assigns the pro¬ 
ton as the missing particle. The intermediate resonances 
to the same pp7r+7r“ final state have been treated as sig¬ 
nals in the tracking efficiency calculation, except those 
with long-lived intermediate states, such as the A parti¬ 
cle. The tight vertex cut used here can result in a lower 
efficiency for long-lived particles, e.g. the tracking effi¬ 
ciency may lower by about 15% at low momentum when 


including the A particle and applying the same selection 
criteria. The total background ratio is less than 0.5% 
in a 3(7 signal region (cr is the missing mass resolution), 
and a similar result can be obtained when assigning the 
missing particle to p, 7r+, or 7r“, respectively. 


etrkjMC) 
^trk (data) 


( 2 ) 


Since n is a subset of iV, the correlation between n 
and N is n, and the covariance matrix of {n,N) will be 
in the form: 


cov{n,n) = {ar,f, cov{N,N) = {aN^ , 
cov{n,N) = {a„y, cov{N,n) = (3) 


with the derivative 

d^^trk ^trk ^^trk 1 / . n 

dN dn '' ^ 

where a„ and CTjv are given by the statistical errors of n 
and iV, respectively. 

The error on the tracking efficiency ) is 


^ ^trk 
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( 5 ) 


In Eq. (ED, due to the independence of etrk{MC) 
and Ctrkidata), the error of the systematic uncertainty 
of tracking ) is 
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HTk(MC) 
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3 Data analysis 

In the selection of the control sample of J/tjj —>• 
pp7r+7r“, only three charged tracks were determined in 
the event selection, leaving the particle under study as 
missed. 


Table 1. Topology analysis when proton is 
assigned as the missing particle using MC sam¬ 
ple of 200M events. 


Channel 

Events 

Ratio 

J/'0 — ?>pp7r+7r“ 

707114 

54.19% 

Resonances to the same final state 

J/'0 — A++ -^ppTT+'K~ 

191565 

14.68% 

J/lp ^ A TT^p ^ ppn'^TT~ 

148190 

11.36% 

J/'Ip ^ A A++ — ^-ppTT+Tr” 

129749 

9.94% 

J/tp ^ PP/q —>■ ppTT+TT” 

33092 

2.54% 

J/ip — >'p7r+A° — >pp7r+7r“ 

29619 

2.27% 

J/ip ^ A°'K~p ^ PP'K'^TT~ 

27464 

2.10% 

J/ip ^ AA —i>p|)7r+7r“ 

13237 

1.01% 

J/Ip ^ pp°p ^ ppTT^n~ 

10524 

0.81% 

J/ip ^ A°A° —i>pp7r+7r“ 

5658 

0.43% 

J /ip^pujp^pp'n^T:~ 

2090 

0.16% 

J/'0 —>■ p /2 (1270)p — ?> pp7r+7r~ 

' 898 

0.07% 

Sum 

1299200 

99.6% 

Other backgrounds 


Sum 

5706 

0.4% 


Charged tracks in BESIII are reconstructed from 
MDC hits. For each charged track, the polar angle 
must satisfy |cos0| < 0.93, and it must pass within 
±10 cm from the interaction point in the beam direc¬ 
tion {\Z\ < 10.0cm), and within ±1 cm of the beam line 
in the plane perpendicular to the beam < 1.0cm), 

which is a vertex requirement. The number of charged 
tracks is required to be at least three. TOF and specific 
energy loss dE/dx of a particle measured in the MDC 
are combined to calculate particle identification (PID) 
probabilities for pion, kaon, and proton hypotheses. For 
the three charged tracks, the particle type yielding the 
largest probability is assigned one by one to each electric 
charge. The invariant mass of one p7r-pair containing no 
missing particle is required to satisfy >1.15 GeV/c^, 
in order to remove A particles from decay chains such as 
J/ij} —>■ AA pp7r^TT~. It is possible some extra tracks 
remain after the selection, and the influence of extra 
tracks is discussed in detail in Sec.lUJv). Figure [T] shows 
the transverse momentum (Pt) distribution of pions and 
protons from the exclusive MC sample after the above 
selection criteria for J/'0 —^-ppTr+Tr”. 


4 Tracking efficiency and its uncertainty 

The number of signal events n or N in Fq. [T] can 
be applied by fitting the missing mass spectrum. The 
A/miss spectrum for proton candidates can be described 
by a Gaussian function and low-degree polynomials, cor¬ 
responding to signal events and background events, sep¬ 
arately. One of the fitting plots from the data 

sample for protons is shown in Figure (Sfa), with the 
requirements Pt £ (0.3,0.35) GeV/c and \cos6\ < 0.7 
(where Pt and cos6 are obtained from missing momen¬ 
tum) . The distribution of Umiss = E^iss - Pmiss (where 
Emiss and Pmisa represent the energy and momentum of 
the missing particle, respectively) is applied for pion can¬ 
didates to achieve the number of signal events, since the 
momentum resolution will cause a discontinuity at the 
zero point on the missing mass distribution of pions. The 
signal and background shape for the Umiss distributions 
of pions are also described by a Gaussian function and 
low-degree polynomials. One of the Umiss fitting plots 
for pions from the data sample is shown in Figure [2](b), 
which requires Pt £ (0.05,0.1) GeV/c and \cos0\ < 0.7. 

The fitting shape used here is inelegant but it is reli¬ 
able, e.g. the signal peak in Figure [2Ka) is not described 
very well by a Gaussian function, but this has little effect 
on the efficiency. The fitting status is similar for both n 
and N, thus their ratio, i.e. the efficiency, is insensitive 
to this minor fitting defect. The number of signal events 
from the exclusive MG sample can be obtained directly 
by counting instead of fitting, and the calculated effi¬ 
ciency is consistent with the fitted efficiency (the differ¬ 
ences are less than 0.1%), which also validates the reli¬ 
ability of the fitting method. Alternative fitting shapes 
have also been tried, e.g. the signal events described by 
a double Gaussian function and the background events 
described by high-degree polynomials, and the differ¬ 
ences in efficiencies are less than 0.1%. Therefore, a 
simpler shape, i.e. a Gaussian function with low-degree 
polynomials, has been chosen for the events description 
to avoid more parameters in the fitting. 

The tracking efficiency and its systematic uncertainty 
are characterized by Pt and cos9, since the efficiency is 
more sensitive to these two variables, which are corre¬ 
lated with the level of track bending and the hit posi¬ 
tions of tracks in the MDC respectively. As an example. 
Figure (SKa) and Figure (SKb) show the one-dimensional 
tracking efficiencies and their systematic uncertainties 
for pions using Pt and cos9, respectively. 

Although different decay processes in physics analy¬ 
ses have different systematic uncertainties for tracking, 
due to the different distributions of Pt and cos9, it is 
possible to provide a general systematic uncertainty for 
various decay processes. If it is represented in two dimen¬ 
sions, Pt versus cos9, then the systematic uncertainty for 
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Fig. 2. (color online) (a) Mmiss fitting status for protons when Pt G (0.3,0.35) GeV/c and \cos9\ < 0.7, and (b) 
Umiss fitting status for pions, when Pt € (0.05,0.1) GeV/c and \cos9\ < 0.7, selected from data samples. The green 
line represents the signal shape described by a Gaussian fnnction, the red line represents the background shape 
described by a low-degree polynomial, and the blue line is their snperposition. 




Fig. 3. (color online) One-dimensional tracking efficiencies and their systematic uncertainties for pions (a) with Pt, 
(b) with cos9 as an example. The red triangles in the upper plots represent the tracking efficiency from the MG 
sample, the black circles in the upper plots represent the tracking efficiency from the data sample, and the circles 
in the lower plot represent the corresponding systematic uncertainty. 


other processes can be obtained by sampling this two- 
dimensional systematic uncertainty. An example of such 
a two-dimensional plot is shown in Figure 0] for protons, 
as a function of cos9 and Pt, in which the colors in each 
box represents the corresponding tracking efficiency or 
systematic uncertainty. 

The results of tracking efficiencies and its system¬ 
atic uncertainties will be helpful for software develop¬ 
ment, especially for MDC track reconstruction. For 
example, one can find the tracking efficiency is less than 
0.9 when Pt < 0.15 GeV/c in Figure |3|a), and the 
systematic uncertainty is also relatively larger. Low 


momentum track reconstruction is difficult, since elec¬ 
tromagnetic multiple scattering, electric field leakage, 
energy loss etc. become stronger, and charged tracks 
with transverse momenta less than 0.12 GeV/c will circle 
within the MDC. The BESIII software group have used 
a TCurlFinder method Q to increase tracking efficiency 
for low-momentum tracks, and improved the MDC MC 
tuning model to decrease its systematic uncertainty M- 
There is still room for improvement of the low momen¬ 
tum track reconstruction and MC tuning, with the assis¬ 
tance of this study. 


010201-4 


















Chinese Physics C Vol. XX, No. X (201X) XXXXXX 



(a) MC (b) Data (c) Uncertainty 

Fig. 4. (color online) Two-dimensional (a) tracking efficiencies from MC sample, (b) tracking efficiencies from data 
sample, and (c) systematic nncertainties for protons as an example. The x-axis is binned in cos6 and y-axis binned 
in Pt, while the color of each box representing a corresponding region efficiency. 




Fig. 5. (color online) Tracking efficiencies from different noise levels from exclusive MC for (a) proton and (b) pion. 
The red triangles represent the efficiencies from a typical low noise level run period (RUN I), and the black circles 
represent the efficiencies from a typical high noise level run period (RUN II). 


5 Validation 

In order to develop a robust methodology for deter¬ 
mining the tracking efficiency, a number of essential vali¬ 
dations are performed as described below, mainly check¬ 
ing the definition, applicability of results, etc. 

(i) Validation of process independent: The tracking 
efficiency should be process independent, which can be 
validated by comparing the tracking efficiencies between 
the single track MC and exclusive MC samples. The sin¬ 
gle track MC is generated by a virtual process for one 
pion or proton track with a uniform phase space, and 
the exclusive MC of the process J/'0 —j-ppTr+Tr" is pro¬ 
duced with a uniform phase space. The momentum and 
polar angle ranges of the single track MC is required to 
be similar to the exclusive MC, then a comparison of the 
tracking efficiency between these two MC samples are 
supplied. The result shows the tracking efficiencies are 
consistent between single track MC and exclusive MC 
samples (the differences are less than 0.1%), which is a 
validation of the process independence. 

(ii) Validation of signal extraction: The inclusive MC 
sample is similar to the data, including both signal and 
background, and the signal is the same as the exclu¬ 
sive MC. The result shows both the tracking efficiencies 
and their uncertainties obtained from the exclusive and 
inclusive MC samples are consistent with each other (the 
differences are less than 0.1%), which proves the signal 


extraction is reliable. 

(hi) Results due to different noise levels: The elec¬ 
tronics noise and beam related background at BESIII 
varies in different run periods during data acquisition, 
so the effects of different noise levels are considered in 
this study. Two run periods were chosen to compare 
the differences in tracking efficiency and uncertainty, one 
being a typical low noise level run period (RUN I), and 
the other a typical high noise level run period (RUN 
II). Figure [5] shows the tracking efficiencies with differ¬ 
ent noise levels from exclusive MC for protons and pions, 
respectively. The results show that high noise level run 
period has about 2% lower efficiency compared with the 
low noise level run period when Pt < 0.3GeV/c. Since 
the different noise levels for different periods of data have 
been introduced into the MC simulation, the system¬ 
atic uncertainties of tracking are not subject to the noise 
level. 

(iv) Comparing results with different vertex cuts: 
Although the qualification for vertex cuts ought to be 
standard, as an elementary requirement which is used 
in almost every analysis, \Rxy\ < 1.0cm, \Z\ < 10.0cm; 
\Rj:y\ < 2.0cm, \Z\ < 10.0cm; and \Rj:y\ < 2.0cm, 
jZj < 20.0cm are three common cuts which can be found 
in different BESIII data analyses. A comparison of these 
three vertex cuts shows that the efficiencies at low Ft 
are obviously influenced by the different requirements. 
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Fig. 6. (color online) Tracking efficiencies from different vertex cuts from exclusive MC for (a) protons and (b) 
pions. The red upward triangles represent the vertex cut \Rxv\ < 1.0cm, \Z\ < 10.0cm; the blue downward 
triangles represent the vertex cut \Rxy\ < 2.0cm, \Z\ < 10.0cm; and the black circles represents the vertex cut 
\Rxy\<2.Gcm, |Z|< 20.0cm. 


Figure [5] shows the tracking efficiencies with different 
vertex cuts from exclusive MC for protons and pions, 
respectively. The results show the efficiency from the 
tightest vertex cut is lower by about 5% compared with 
the loosest cut at low Pt- The data sample behaves in 
the same way as the MC sample, so the systematic uncer¬ 
tainties of tracking efficiency are insensitive to the vertex 
cuts. 

(v) The influence of extra tracks: The n and N in 
Eq.[T] were explained to be the number of signal events 
in which the number of reconstructed charged tracks are 
greater than three and four, respectively, which allows 
some extra tracks in the definition. A reconstructed 
track with some extra fake reconstructed tracks was 
introduced to the tracking efficiency in this study, but 
the extra tracks may be avoided in most physics analyses, 
on account of the requirement of zero net charge. The 
influence of the extra tracks should be evaluated. When 
count the number of charged tracks in each event, only 
a few events have the number of charged tracks greater 
than four, which means the extra fake tracks are seldom 


reconstructed. The tracking efficiencies with or without 
the extra tracks were compared for both MC and data 
samples, and found to have consistent tracking efficien¬ 
cies (the differences are less than 0.1%). Therefore, the 
influence of extra tracks can be neglected. 

6 Summary 

A robust methodology for determining the tracking 
efficiency and its systematic uncertainty is presented in 
this paper. A clear definition of tracking efficiency and 
its systematic uncertainty was given, a clean sample of 
J/Ip ^ pp7r+7r“ was selected and a number of essen¬ 
tial validations performed to guarantee the reliability of 
the method. The tracking efficiency and its systematic 
uncertainty was determined in both one and two dimen¬ 
sions as a function of both transverse moment and polar 
angle, which were shown as examples. It is possible to 
use the general results of tracking efficiency and its sys¬ 
tematic uncertainty for the analysis of other processes. 
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